clear all
set more off
set maxvar 32000
set matsize 11000
*************************************************************************************************
*set scheme aspen, perm
grstyle init
grstyle set plain, horizontal grid

use "$folder/data/data_ready.dta", clear
cd "$folder/output/tables"
local x=18

xtset state year, yearly delta(2)

label variable growth_gsp "Growth GSP"
label variable log_real1gsppc "Log Real Income per capita"
label variable growth_pc "Growth Real Income per capita"
label variable lestab_D "Growth Number Estabilishments"
label variable lprofits_mill_D "Growth State Profits, Millions"
label variable lcbp_emp_D "Growth Employment"

label variable log_provisions_18 "Legislative Detail"
label variable log_cond_list_prov_18 "Contingent Provisions" 
label variable log_non_cond_list_prov_18 "Non-Contingent Provisions" 
label variable shock_prov18 "Shock to Legislative Detail"
label variable shock_cond_list_prov18 "Shock to Contingencies"
label variable shock_non_cond_list_prov18 "Shock to Non-Contingencies"

label variable gsp_initial "Initial Gross State Product"
label variable growth_initial "Initial GSP Growth"
label variable gsp_pc_initial "Initial GSP per capita"
label variable lestab_D "Growth Number Establishments"
label variable lprofits_mill_D "Growth State Profits, Millions"
label variable lcbp_emp_D "Growth Employment"

*******************
* FS-OLS-RF
*******************

* First Stage
reghdfe log_provisions_`x' shock_prov`x' , absorb(state year) cluster(state)
outreg2 using Table_2_FS-OLS-RF, excel replace label keep (shock_prov`x') ctitle(First Stage)   nocons addtext(State FE, X, Time FE, X)  alpha(0.01, 0.05, 0.1) symbol(**, *, +) nor2  nonotes addnote(SE clustered by `e(clustvar)', **p<.01; *p<.05; +p<.1) addstat( R-squared,  e(r2)) 
* First Stage + state year trends
reghdfe log_provisions_`x' shock_prov`x'  , absorb(state##c.year year) cluster(state)
outreg2 using Table_2_FS-OLS-RF, excel append label  keep(shock_prov`x')  ctitle(First Stage)    nocons addtext(State FE, X, Time FE, X, State-Specific Trends, X)  alpha(0.01, 0.05, 0.1) symbol(**, *, +)  nor2 nonotes addnote(SE clustered by `e(clustvar)', **p<.01; *p<.05; +p<.1) addstat( R-squared,  e(r2))

* OLS
reghdfe growth_pc log_provisions_`x' , absorb(state year) cluster(state)
outreg2 using Table_2_FS-OLS-RF, excel append label keep (log_provisions_`x') ctitle(OLS)    nocons addtext(State FE, X, Time FE, X) alpha(0.01, 0.05, 0.1) symbol(**, *, +) nor2  nonotes addnote(SE clustered by `e(clustvar)', **p<.01; *p<.05; +p<.1)  addstat( R-squared,  e(r2))
* OLS + state year trends
reghdfe growth_pc log_provisions_`x' , absorb(state##c.year year) cluster(state)
outreg2 using Table_2_FS-OLS-RF, excel append label  keep(log_provisions_`x')  ctitle(OLS)    nocons addtext(State FE, X, Time FE, X, State-Specific Trends, X)  alpha(0.01, 0.05, 0.1) symbol(**, *, +)  nor2 nonotes addnote(SE clustered by `e(clustvar)', **p<.01; *p<.05; +p<.1) addstat( R-squared,  e(r2))

* Reduced Form
reghdfe growth_pc shock_prov`x' , absorb(state year) cluster(state)
outreg2 using Table_2_FS-OLS-RF, excel append label  keep(shock_prov`x')  ctitle(Reduced Form)    nocons addtext(State FE, X, Time FE, X) alpha(0.01, 0.05, 0.1) symbol(**, *, +)  nor2 nonotes addnote(SE clustered by `e(clustvar)', **p<.01; *p<.05; +p<.1) addstat( R-squared,  e(r2)) 
* Reduced Form + state year trends
reghdfe growth_pc shock_prov`x' , absorb(state##c.year year) cluster(state)
outreg2 using Table_2_FS-OLS-RF, excel append label  keep(shock_prov`x')  ctitle(Reduced Form)    nocons addtext(State FE, X, Time FE, X, State-Specific Trends, X)  alpha(0.01, 0.05, 0.1) symbol(**, *, +)  nor2 nonotes addnote(SE clustered by `e(clustvar)', **p<.01; *p<.05; +p<.1) addstat( R-squared,  e(r2))

